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ABSTRACT 


A full-scale model of a double-wedge, simulating a jet 
vane used in guided missile thrust vector control (TVC), was 
constructed for subsonic wind tunnel testing to evaluate the 
feasibility of utilizing a thermal imaging (infrared) system 
vice thermocouples to determine a surface heat transfer coeffi- 
cient. The model was preheated, then allowed to cool via 
forced convection to ambient conditions. Surface temperature 
readings were taken using thermocouples and a thermal imaging 
system during the cooldown process. These readings were then 
reduced to solve for a heat transfer coefficient (h), using 
Newton's Law of Cooling. The results were compared with the 
theoretical solution for a flat plate and with the results 
obtained by testing an actual TVC vane configuration under 


mm ident.cal flow conditions. 
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LI.  INTROBUETTON 


A. BACKGROUND 

The use of thrust vector control (TVC) in future cruise 
missile guidance and control systems is presently being inves- 
tigated. Removing external aerodynamic steering surfaces 
(fins) would allow installation of a single weapon configuration 
in multiple platforms with little or no restructuring of each 
platform's unique launch capabilities, thereby reducing desi $ 
development, and construction costs. In addition, the maneu- 
vering envelopes required of present and future missiles may 
contain regions in which external aerodynamic surfaces are 
inadequaces fOrscon Erol pot 

The incorporation of TVC systems for missile steering re- 
quiresithe design of jet vancs Or Ot mere I lc e vi D 
that are mounted in the rocket exhaust flow. Such devices can 
be rotated to varıous angles of attack in order to redirect the 
flow, thus maneuvering the missile in the desired direction. 
The vanes must be made to withstand severe steady and transi- 
ent thermal loads, shock wave impingement, and erosion from 
motor exhaust particulates during extended use. 

One of the first steps in the design of such a vane is to 
develop a model of the thermal environment to which the vane 
will be subjected. This includes the determination of surface 
heat transfer coefficients for particular vane configurations 


under forced convection Conditions. 


The following discussion concerns the development of 
experimental methods to determine a heat transfer coefficient 
on a heated surface in subsonic flow. Utilizing state-of-the- 
art infrared technology, in place of thermocouples, to determine 
surface temperature gradients as a body is heated (or cooled, 
as in this case) would reduce the time and effort required to 
Sinstrüm m e models, The method also has the potential 
to reduce or eliminate the need for elaborate data acquisition 
systems, while maintaining an accuracy that is at least as 
good as that obtained from conventional methods. Finally, 
the methods described here provide a global as well as local 
view of surface heat transfer. Thus it is possible to detect 
and analyze regions of particular interest that might otherwise 
go unnoticed in the data produced by point temperature 


measurements. 


BP OBJECTIVES 
The objectives of this study are to: 


1) Design a typical: TVC model instrumented for surface 
temperature measurement using both thermocouples 
and thermal imaging technology. 


2) Record surface temperatures of the model in a subsonic 
wind tunnel environment. 


3) Compare the results from both measurement techniques in 
calculating the heat transfer coefficient, using 
rr Sx law of Cooling, and correlating the results 
with the theoretical solution for a flat plate under 
me lí j cal Flow conditions. 


4) Utilize the thermal imaging technigue to analyze the 
surface temperature map of an actual TVC vane 
configuration. 


5) Present conclusions and recommendations for future 
aplica ions in this area. 


II. MODEL DESIGN AND CONSTRUCTION 


A. MODEL DESIGN 

A double-wedge design was chosen for this experiment for 
several reasons, including ease of construction, similarity 
to typical TVC configurations used in previous research, and 
sufficient similarity to a flat plate to allow interpretation 
of results. A drawing of the wedge is included as Figure l. 

The dimensions of the model were chosen to reflect essen- 
tially full-scale, but constrained by the necessity of having 
a Biot number less than 0.1. This constraint is mandated by 
the use of the lumped-heat capacitance analogy in this experi- 
ment to determine the heat transfer coefficient (Ref. 11. 
This method of analysis assumes that the internal resistance 
of the body is negligible compared to the external resistance 
and, therefore, a uniform temperature distribution exists 
throughout the model. By ensuring that the model has a Biot 
number less than 0.1, this analytical method will produce 
reliable results. 

The convection heat loss from the heated model is assumed 
to equal the decrease in its internal energy. Thus, the 


energy balance expressing the rate of heat transfer, q, is: 


where: 


10 


h = heat transfer coefficient, BTU/hr-ft^-ep 


surface area, E 


ne - 
T = model temperature at time t, °F 
e ww ce ercan temperature, °F 
c = heat capacity of model material, BTU/lbm-°F 
V = model volume, ec? 
EE ere En 1. BTU/hr-ft&-?F 
R s 3 
oa. MEEN POU EE. 
Seen the initial condition of T = To dt time t = 0 the 


solution to the differential eguation becomes: 


DEP 
T 7 = exp ~- (Bi xFo) 
O CO 
where: 
I | 3 

Bi Biot number = KA 
T EA. 2 

FO Fourier number = sc y. Ë 


Dr CONSPRUCTION 

Copper was chosen as the model material due to its excel- 
lent thermal conductivity (k = 230 BTU/hr-ft-°F) and reasonably 
good machinability. With the dimensions as shown in Figure 1l, 
giving V/A = 0.01297, the surface heat transfer coefficient 


2 
could theoretically rise to well over h = 1400 BTU/hr-ft -°F 


mr 


before the restriction on Biot number would be violated. As 
the expected heat transfer coefficient was less than 60, these 
design decisions ensured excellent accuracy in applying the 
method. 

The copper wedge was instrumented with copper-constantin 
(type T) thermocouples, 30 gauge. The thermocouples were 
placed as shown in Figure 1, with one of the thermocouples 
embedded at depth of 0.30 inches to detect any distinguishable 
temperature difference between the surface and the interior. 

As can be seen in the subsequent tables and graphs, the assump- 
tion of Newtonian cooling proved to be accurate. The wedge 

was also covered with a thinzeasr Ot flat black Tait co ME 
.crease the emissivity of the surface to a minimum of 0.95 

[Ref. 1] in order to enhance the thermal image detection of the 
infrared systen. 

The thermocouples were attached to a dual-channel strip 
chart recorder, and tota tata len peras e IE 
digital readout was used to obtain reference values of wedge 
temperature during heatup, and free stream temperatures during 
testing. All instruments and thermocouples were calibrated 
from 32°F to 150°F using the Department of Mechanical Engineer- 
ing calibratron facility. 

The wedge was then mounted on two aluminum stanchions, 
attached to the plexiglass insert in the floor of the wind 
tunnel test section. Within the plexiglass, a six-inch diam- 


eter hole was cut for insertion of the infrared-transparent 


12 


window through which the thermal imaging system could view 

the wedge surface. This arrangement is shown in Figure 2. 
Plexiglass sidewalls were added to form a thermal barrier to 
attempt to negate any spanwise heat transfer from the wedge 
sides, and to form a channel for uniform air flow past the 
wedge. The model was heated by using an electric pad, capable 
of raising temperatures from ambient to approximately 140°F 


within a thirty minute time frame. 


RES 


III. THERMAL IMAGING SYSTEM 


A. DESCRIPTION 

The thermal imaging system used in this procedure is a 
Probeye Thermal Video System, e es 300, ane rcc s 
Hughes Aircraft Company [Ref. 2). This system provides real- 
time temperature maps of objects through 16 color bands dis- 
played on a color monitor. The temperature range of the system 
is from -4?F to 536°F edge eo 230°C). Bach m ie 
scale can represent temperature sensitivities from 0.5 to 
30 HU SE 20°C), depending upon the desired overall tempera- 
ture range. Emissivity adjustments can be made from 0.0 to 
1.0 in increments of 0.01, in the spectral wavelength band of 
2.0 to 5.6 umeters. The unit includes a movable cursor which 
can be placed at any location on the viewed surface, giving a 
digital readout of the temperature at that location, and a 
real time clock/calendar precise to 0.1 seconds for determining 
time-variant conditions. 

The use of this imaging system also required the use of 
an infrared optic window for viewing the model while maintain- 
ing tunnel integrity. A six-inch diameter, two millimeter 
thick window fabricated from magnesium fluoride was obtained 
from Eastman Kodak Company. This type of material, tradename 
IRTRAN 1, was chosen to match the optical wavelength charac- 


teristics of the imaging system. As manufactured, the window 
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allowed over 95% transmittance in the 2.0 to 5.8 umeter 
range, while capable of withstanding pressure differentials 
oil atmosphere at 25°C, withwa factor of safety of four 


(Ref. 3]. 
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IV. EXPERIMENTAL PROCEDURE 


A. WIND TUNNEL TESTING 

The model assembly, as shown in Figure 2, was installed 
in the test section of the subsonic wind tunnel. This wind 
tunnel is a pull-down system capable of producing test section 
free stream velocities of up to 335 feet/second by means of a 
variable speed fan located at the tunnel exit. The free-stream 
turbulence level of the tunnel (0.2%-0.5%) was deemed suita- 
ble for the purposes of this experiment [Ref. 4]. 

At the start of each test, with the tunnel fan off, the 
wedge was wrapped with the heating pad, and brought up to 
temperatures of approximately 130°F. The pad was then removed, 
wind tunnel access doors closed, and the fan started. Runs 
were made at 50% and 75% fan speed, corresponding to 200 and 
306 feet/second, respectively. Free stream velocities were 
calculated using the installed manometer, calibrated in inches 
HO. Free stream temperatures were determined by both a 
calibrated thermometer and a thermocouple inserted in the flow. 

The cooling process was then recorded using the installed 
thermocouple output to the strip chart recorders, and the ther- 
mal imaging system connected to a video tape recorder for 
simultaneous data acquisition. In initial tests, the tempera- 
ture sensitivity settings on the imaging system were selected 


at 5.0°F, which yielded an 80 degree range from 50°F to 130°F, 


dié 


Ee view the entire cooling process. In subsequent tests, a 
sensitivity of 2.0°F was selected, giving a 32 degree range 
from 90°F to 122°F. This lower sensitivity/range setting was 
necessary to obtain a good agreement with the thermocouple 
data. 

Serene Ner able amount of system testing, this process 
was undertaken four times at each of the two fan speeds for 


data evaluation and correlation. 


ELM DATA REDUCTION 

Applying the formula for the lumped-capacitance analysis, 
the transient data from the thermocouple readings and the 
thermal imaging system were used to deduce the corresponding 
heat transfer coefficients. The calculations were performed 
using TK! SOLVER, a microcomputer spreadsheet [Ref. 5]. Tables 
l and 2 exhibit the program equations and variables, which 
provide not only the heat transfer coefficient, but also an 
estimate of the uncertainty of that calculation, and the 
corresponding Biot number, for each of the two sensor types. 
The uncertainty analysis of the heat transfer coefficient 
calculation is included as Appendix A. 

The reduction of the thermocouple data was fairly straight- 
forward, as the output was temperature versus time. However, 
reducing the thermal imaging system data required running the 
video tape recording frame by frame (each frame corresponding 


to 1/60 second) to determine the temperature at a given time. 


di 


This method generated the largest uncertainty in the calculation 
because the visual color bands and digital readout of the sys— 
tem changed only when the temperature of the body changed 

more than the sensitivity setting. For example, at a 5.0°F 
setting, the readout would persist for perhaps three ccm 

of test time before changing from LOSE to. 105°F fh legen 
reality, as displayed by the thermocouple readings, the surface 
temperature could be any temperature within this range ata 
particular instant. (This problem is sous to that due 

to a large interval size in a finite-difference calculation.) 
There also existed occasions when, particularly at temperatures 
near thermal eguilibrium, the readout would take relatively 
long periods of time to stabilize, which made determining the 
exact temperature/time data difficult. However, these diffi- 
culties were overcome by operator experience and use of the 
narrower sensitivity setting for correlating data in smaller 


segments rather than over the entire cooling span. 
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NANA ERIN IDA RESULTS 


A. FLAT PLATE CORRELATION 

The analysis of a flat plate of identical material, length, 
and width in identical flow conditions was used to attain a 
degree of confidence in the accuracy of the results of the 
wedge experiment. This analysis is based upon determining 
the Reynolds and Prandtl numbers for the flow, and solving the 


ww Da equation Hor the average, Nusselt number: [Ref. 6] 


1775 > 


Nu — 0.664 Re Pr 


mie Nusselt number was then used to determine the average 
IE MA Ster coerrıcıent from: 


Nu k 
H 





SE 
N 


where 
VEEL OF flat plate, ft. 


These procedures resulted in a heat transfer coefficient of 
22 Emi me op for free stream velocity of 200 ft/sec, and 
28 BTU/hr-ft--°F a eros oF 3060 ft/sec. The results of 


this calculation are shown in Appendix B. 


ES 


B. THERMOCOUPLE RESULTS 

Using the strip chart recorder traces, the heat transfer 
coefficient at the surface thermocouple location was determined 
from the lumped-capacitance solution of exponential temperature 
decrease with time. For the free stream velocity of 200 ft/ 
sec, h varied from 332 BTU/hr-ft^-*F, with an average 
uncertainty band of Ł1.75 BTU/hr-ft^-°F. At vetoc seyi 
306 ft/sec, h had values ranging from 41.03 to 42.7 BTU/ 
E with average uncertainty of +1.80 BTU/hr-ft“-?r. 


These calculations are shown in Tables 4-ll. 


C. THERMAL IMAGING SYSTEM RESULTS 

As shown in Tables 12-19, using the same analysis and 
temperature range, h values were determined to be 27.7 to 
3585 BTU/hr-ft^-*F, with an band of £1229 BYU; 
hr e oe for the lower free stream velocity, while values of 
39.9 to 41.8 were obtained at higher speed. Figures 3 and 4 
show the comparison of wedge surface temperature versus time 
for the two sensors at each velocity, and Figures 5 and 6 
show e ed values from each sensor versus time, again 
at each speed. 

It is important to note that this measurement corresponds 
to conditlons at a particular point on the wedge, that is, the 
location of the corresponding thermocouple. Further testing, 
with the cursor positioned at other points, would allow the 
development of a complete heat transfer map of the wedge sur- 


face. However, viewing of the entire wedge surface showed that 
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the transient isotherm variation was essentially uniform over 

a large portion of the surface (with the exception of narrow 
bands in the three-dimensional boundary layer regions near the 
sidewalls). Figure 7 is a still photograph of the wedge taken 
from the video output, which clearly exhibits this condition. 
Thus it may be deduced that the values determined above, at 

the thermocouple.location, are representative of conditions 
over a major portion of the wedge surface. The ability to make 
such observations is a significant virtue of the infrared 


thermography method. 


D. TVC VANE ANALYSIS 

An actual TVC vane. was received from NWC China Lake, Ca 
for analyzing the effect of complex geometry on convective heat 
transfer using the thermal imaging system. 

This vane, shown in Figure 8. was constructed of copper- 
impregnated tungsten (803 W, 203 Cu), with the following 


material properties: 


lbm BTU BIU 


D — 1076.544 s k = 148 EES EER C — 0:037 NM "vr 
ft 
=> 3 
The vane volume was calculated to be 2.295 x10 ne wat) 
- 2 
a corresponding surface area of 1.5551 x10 i ft . All values 


were then entered into a formula sheet as before, shown in 


Table 3. 
The vane was mounted on two aluminum stanchions with the 


same baseplate used by China Lake in their experiments. The 


SE 


thermal imaging system was set up to view the vane through the 
Side port of the subsonic wind tunnel with the IRTRAN optic 
window installed. The identical heating procedure was used, 
however, for reasons unknown, the pad would only reach 105°F. 
This was deemed sufficient, as it still gave a 40 degree tem- 
perature differential. The wind tunnel was operated at 50% 
flow speed (200 ft/sec). 

Three test runs were performed, with the cursor at the 
midspan for runs l and 3, and on the leading edge at the base- 
plate boundary for run 2. The cursor was moved for run 2 to 
investigate and analyze a more rapid temperature gradient 
noticed in that area in the previous run. Run 1 yielded an 
average h of 41.5 BTU/hr-ft?-or +3.23, run 2 yielded an average 
Of DA and “ri gave an average of 42.4 38.9, The 
largest Biot number was calculated as 0.0063, well within the 
lumped-capacitance analysis criteria. 

Figure 9 is a color thermograph of the vane during cool- 
down. As compared to the thermogram of the double wedge 
(Figure 7), the temperature map is fairly uniform over the 
body, with the exception of faster cooling at the leading and 
trailing edges, and an interesting area near the top of the 
vane below the apex. While vane thickness would explain the 
shape of most of the isotherms, the upper region undoubtedly 
is affected by a unique flow field, which increases the heat 
transfer coefficient by approximately 83 (Table 22). Tables 
20-22 show the results of these tests, and Figures 10 and ll 


are graphs of temperature and h versus time, respectively. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


AE ONCLUSIONS 

The results obtained from the thermal imaging system are, 
within a small experimental uncertainty, in agreement with 
the results of the thermocouples for all conditions tested 
here. Both experimental method's results were also suffi- 
ciently close to the theoretical flat plate solution to allow 
confidence in their accuracy, primarily due to the small wedge 
angle used. 

The primary disadvantage cf the thermal imaging system is 
the increase.in size of the uncertainty band as the tempera- 
ture sensitivity is widened. For large sensitivities of 5 
degrees or more, Ss results become comparatively less accurate, 
again similar to using large step sizes in finite element 
calculations. Future users of this system in experiments 
similar to this one will have to adjust the sensitivity settings 
with respect to the overall temperature range desired, most 
likely by trial and error, to get the degree of accuracy 
reguired. 

Nevertheless, the results obtained in this experiment 
lead to the conclusion that the thermal imaging system, when 
properly employed, would be a very good alternative to thermo- 
ou Is applications similar to this one. It provides 


accurate local and global readings with no need for 
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complicated data acguisition/reduction eguipment and extensive 
instrumentation. The system is transportable, very reliable, 
and easy to operate. Although the system reguires special 
infrared optic material for viewing models in an enclosed 
environment, several vendors were found who are capable of 
providing the material in various shapes and sizes at a reason- 


able expense in time and money. 


B. RECOMMENDATIONS 
The findings of this research lead to a recommendation for 
using this system for surface temperature measurements in the 
next phase of this project, using the wedge model in a super- 
sonic wind tunnel for the same heat transfer coefficient calcu- 
lation. This model was designed for ease of installation in 
the supersonic tunnel, when it becomes available, and an optic 
window has been obtained that will allow infrared viewing of 
the model while withstanding the environment (temperature and 
pressure differentials) in supersonic flow. Particular problems 
that might be encountered in supersonic applications include: 
1) Sufficient difference in model and free stream 
temperatures to cause adequate heat transfer--the 
model may have to be heated prior to testing by 
some internal means (installation of wires for 
Joulean heating) to increase the recovery factor. 
2) The on-site supersonic tunnel has a significantly 
smaller run-time available due to its pressure 
vessel configuration, thus yielding short run 
durations (a few minutes) and long (hours vice 
minutes) pressurization durations between tests. 
3) As other model shapes and materials are considered, 


the necessity for maintaining the Biot number less 
than 0.1 could become critical, or another method 
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of analysis would have to be chosen if the lumped- 
capacitance assumption is not valid. 


It is also recommended that the next experimenter familiar- 
ize himself with the system thoroughly, and invest time in 
viewing the video tape outputs of the system to get experience 
Bere mperatcure determination. This experience will prove to 


be a great factor in enhancing the accuracy of the results. 
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APPENDIX A 


UNCERTAINTY ANALYSIS CALCULATIONS 


= e HAL : E 
Oo = T g = (sg GeV. Formula used in caleuslat mer 
h = - Anðlpc V) Solvumndgesor h 
At 
h — B(O,t ALME) 
2 oh 2 oh 2 oh 2 oh 2 
Ah = EET AO) Gei (SE At) + (7 AA) + (37 AV) 
E72 
oh 2 oh 2 
t (55 Ap)” + (57 &c)"] 


where A's are the uncertainty of the measured quantities. 


mw 
oh — DOM RR; 2 E 2 
(1) 36 AQ = SAT. | SC AD) EE (5 AN) ] 
where: 
2 2 1⁄2 
Nis Tans Dam T AM p r u) ; 
2 2 1 2 
AD. = AT 0 ER 
AT = 2.0 °F for IR data; "AT = OR Fore NEO OUR 


data 


26 


(2) ən e e 
At 


where: 
ACO seconds = 2.7778 E -5 hours 


from IR system timer precision. 


an ma O pe V 


(3) — AA = AA 
JA Act 
where: 
ee p 
AA = 00000 in = 6.9444E-6ft 
from machining precision. 
3h "INE H 
(4) ay AV = e AV, 
where: 
Ce 3 
AV si EA OES in —-. 5.7734 1) ft 
mFronmemacmınıng precision. 
oh won e V 
(5) E Ap = e a âp, 
where: 
Ap = 6.9082 U (ref. 71 
ft 
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oh sn EE 


(6) De AC M — TI SCORES Ac, 
where: 
B BTU 
AC — 0 001 Tbm-°F [Ref. 7] 


All above equations were included in the data reduction 


programs to solve for the uncertainty in h (delh). 
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APPEIIEIET B 


FLAT PLATE SOLUTION 
Given: T = 60°F = 520°R 1"H.0 03 I EDS 


ph rom the empirical flat plate in parallel flow 








equation: 
e occa Re; / Í Dr 
freer. 6] 
— o Nuk 
Ds D 
(JE) Re, : Nor — 520°R: (Ref. 7] 
_ lbm E = BTU 
Der = 12 205395=5 EE EE Pr = DEN k = 0.01561 e 
EE 0.0743 m e 
alr 3 
Le 
32 E DM (8.9 H40) (0.0361) (64.348) (144) 
D Y - 0.0743 
2 ft^ ft | 
va. 40080 4 — UNE. 200:2 —- “for experimental 
S MESI speed flow 
EE : 
Poren = 20.8” HO, oa 0o s for experimental 


high speed flow 
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YU x (007 2) (200.2 TO. 2063) 





Re. = BEE = ==) 000 — ISE 
(0.0743) (506) (0. 20559 
Re. OE eme ae EAR 
Assuming Re s = 5.0E5, laminar flow conditions exist in ka 
cases. 
(2) Nu: For low speed flow: 
Na = 0:664(2.574E5) ^ 2008970) ^ ^) = 300.5 
For high speed flow: 
en / 
Nu EE 
(3) E For low speed flow: 
¡A GE e 
: hr-ft -?F 
For high speed flow: 
p. Nik _ (371.5)(0.01561) , 7. BTU 
P i hr-ft^-?F 
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FORMULA SHEET FOR THERMOCOUPLE 


id input 
L 
L 
L 
L 
7771717178 
n 
550 
20941 
.00346013 
L 
250 
L 
L 


Name 


theta 


TE 
RE 


delh 


a 


TO 


APPENDIX C 


TABLES 


TABLE 1 


output [Unsere 


DEN 


ai ep 


Deu 


hr-ft2-°F 


res 2 
hours 
ibm, Ee 3 
btu/lb-F 
Tias 
Its 


DE 
hr-ft-°F 


Seconds 


Du 


DATA REDUCTION 


Comment 

dimensionless 
temperature 
temperature at time t 


free stream 
eemperchelir e 


initial wedge 
temperature 


heat transfer 
coefficient 


uncertainty in h value 


wedge wetted perimeter 


copper density 
specific heat capacity 
wedge volume 

heat transfer 


wedge conductivity 


(dh/dtheta)* (delta theta) 
çimde) (delta) 
(dh/dA)* (delta a) 
(dh/dv) * (delta V) 

(an dro) * (delta ro) 

ar de IE (delta e) 


TABLE 1 (CONTINUED) 


S Rule 
* theta = Ben 
ochebra DINE c a 
a = a en) 
> Brom. Ten 
* time = t*3600 


*u —İre e V) theta*at) (TE Tinf |) Ain oD 
(1.655697 Wr mamo) 2) 


* VV = 2E 7116-5 9 m mta) ro*c*v/ att c2) 
*w = (ln(theta)*ro*c*v*6.9444e-6/(a^2*t)) 

* x sd eta) “Fo AE je-13/ (a*t)) 

* y NEC Eheta)*c*VUu RAAI E)) 

* z = (lin(theta)*v£ro> wer (a>e)) 

* delh = sort (u”"2+vv"2+w"2+x"2+y"2+z"2) 
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ABLE 2 


FORMULA SHEET FOR INFRARED SYSTEM DATA REDUCTION 


Sin put Name output Une Comment 
L theta dimensionless temperature 
L EE E temperature at time t 
ELE F ‘free stream temperature 
PE F initial wedge temperature 
L h | — btu heat transfer coeffi- 
hr-ft?-op cient | 
L delh GN üneercainty in yele 
717718 a BEAST m wedge wetted perimeter 
L js hours 
558 ro lom Apes copper density 
pal G |o Eu E specific heat capacity 
.00346013 V TE wedge volume 
L g Dee ne heat transfer 
290 k o btu = wedge conductivity 
Rer- re EF 
L BTOT 
L time seconds 
u (dh/dtheta)* (delta theta) 
vv de)” (delta t) 


Ww (dh/dA)* (delta a) 
x any (delta) 
y dido) “(delta Ho) 
Z (An de) * (delta c) 


S Rule 
wire Timi) 
WEE ET 
= Kat 
theta = exp (Tr) 


* q = h*a* (Tt-Tinf) 


> 


TABLE 2 (CONTINUED) 


hv 


time - t*3600 


u s(ro*c*Vilitheta EEGENEN 


+ CLOS MU DIS p 5S RN 


vv = (2.777777e-5*In(theta) *ro*e*Vv/(a*t”2)) 


w = (ln(theta)*ro*c*v*6.9444e-6/(a^2*t)) 
x = (1n(theta)*ro*c*5.787e-13/(a*t)) 
y =wüm(Lheta)íxc+*V*6,.05° E)) 


N 
II 


(ini(theta) *V*ro <a a) ) 


delh = sgrt(u”2+vv"2+w"2+x"2+y"2+z"2) 
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TABLE 3 


FORMULA SHEET FOR TVC VANE DATA REDUCTION 


Rule 
(Tt-Tinf) 
* AA 
SC AE] 
-(h*a*t) 


theta = exp( ) 


(ro*c*V) 


e Ha (le bime) 


mm nu 
time - t*3600 


Mere OY (theta*a*t)) *((L.4142*(Tt-Tinf)/(Ti-Tinf) 2) 


we (2%, 2360 ('Pi-Tinf))^2) 


vv = (2.777777e-5*1n (theta) *ro*c*V/ (a*t*2) ) 


w = (ln(theta)*ro*c*V*6.9444e-6/(a"2*t)) 
ro (in (Eheta) *ro*e*5 .787e-l13/ (a*t)) 

NE nUortheta)*c*v*6.9082/(a*t)) 

Ad In theta) “V“ro*,.001/(a*t)) 

MEM sgrk (un 2T-Tvv 2tw 2+x 2ty 2tz 2) 
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TABLE 3 (CONTINUED) 


St input Name Output Umi Comment 
L theta dimensionless temperature 
L oI Die F temperature at time t 
61 Tore F free stream temperature 
97 Tr F initial TVC vane 
cemperacure 
Di E 
L h —— heat transter coerrie ms 
hr-ft^-»p 
L delh uncertainty min Mva ME 
155507 a fe 2 TVC vane surface area 
L e hours 
1076.544 ro lbm/ft^3 TVC vane density-80$ W, 
20$ Cu 
07 C be LEER specific heat capacity 
UA Vv fees TVC vane volume 
L g btu/hr heat transferred 
Sea T 
148 | k ir SOP TVC vane conductivity 
L BIOT 
L 4753 time seconds 
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TABLE 4 


RESULTS OF THERMOCOUPLE CALCULATIONS, LOW TUNNEL SPEED 


kum: nl 
Duel Air Velocity: 201.3 -Ít. 
Sec 
Average Biot Number: 0.00164 
Average Uncertainty in h: m 20 — 
hr-ft -°F 
Time (Seconds) Temperature (“F) Ba EE ) 
= ——— LAA hr-ft-°F 
050 el 725 
6.7 112-5 30.92 
14.5 10725 30200 
2:27 1025 BIG 
220 97:5 30.42 
43.0 925 o IE 
5520 867.5 SE] 
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TABLE 5 


RESULTS OF THERMOCOUPLE CALCULATIONS, LOW TUNNEL SPEED 


Run: #2 
| «^ Fe 
Tunnel Air Velocity: 193.3 —— 
sec 
Average Biot Number: 000156 
PEU 
Average Uncertainty in h: 214,95 5 
| r= r 2m 
Time (Seconds) Temperature (°F) hår peu ) 
Dr te} 
O 0 11455 
Eu 112 5 24. 32 
5.6 11095 30299 
Ir 10825 28-56 
125 TOS 20 324 
16.0 104.5 28:85 
21006 10205 26:97 
24.6 100m5 27.48 
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TABLE 6 


RESULTS OF THERMOCOUPLE CALCULATION, LOW TUNNEL SPEED 


Run: #3 
¡Manel Air Velocity: u s 
sec 
Average Biot Number: 0200159 
; Z Deu 
Average Uncertainty in h: +y 1555 5 
Dr tas FE 
Time (Seconds) Temperature (°F) h ( SS ) 
— A AAA ——— u Mr TE Fr 
0.0 AO 
00 | PE ; 28.70 
20.0 07,0 29.47 
200 1017 EA 
40.0 96.7 29.92 
50.0 92.8 29.54 
60.0 89.2 29954 
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RESULTS OF THERMOCOUPLE 


TABLE 7 


CALCULATIONS, LOW TUNNEL SPEED 


Run: #4 
Tunnel Air Velocity: 2 ss" 
sec 
Average Biot Number: 0.00182 
Average Uncertainty in h: zl — 
hr-ft -°F 
Time (Seconds) Temperature (°F) hut Se ) 
=> — — hr-ft-*F 
0.0 11190 
5.0 3 0p 32.44 
10.0 104.3 34.23 
1.5.30 101.4 0 12 
20.0 BT Ba] 
2550 96.2 33469 
3020 557 3992 
3520 91.6 gio 
DEN Qu 2 
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TABLE 8 


RESULTS OF THERMOCOUPLE CALCULATIONS, HIGH TUNNEL SPEED 


Run: #5 
Tunnel Air Velocity: 308.2 tt 
sec 
Average Biot Number: 00022 
DEL 
Average Uncertainty in h: t 17838 5 
hear ri 
Time (Seconds) Temperature (°F) h(—Ptu y 
— ——— A hr-ft-°F 
0.0 119.0 
Stee 114.3 40.37 
TOO 10927 41.45 
MO 1:02:85 41.73 
20.0 101.7 41.95 
25.0 98.3 42.03 
30.0 Í 42.00 
DL 92272 42.26 
5 50 Ou) 42.66 
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TABLE 9 


RESULTS OF THERMOCOUPLE CALCULATIONS, HIGH TUNNEL SPEED 


Run: #6 
Tunnel Air Velocity: 30 s 
sec 
Average Biot Number: 0.0022 
Average Uncertainty RM. Zİ - 
hao SSE 
Time (Seconds) Temperature (“F) neun N 
Ë — = wo å ooo EE EE 
0.0 11550 
2.0 Lis | 38 15 
TOTO 712003 38.64 
15256 1067 37.48 
20) DOE 40.84 
25.0 3555 AZ 
30.0 oi 41.07 
25.0 92.4 41.38 
900 ENG 41.54 
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TABLE 10 


RESULTS OF THERMOCOUPLE CALCULATIONS, HIGH TUNNEL SPEED 


Bun: #7 
Tunnel Air Velocity: SNO Lt 
sec 
Average Biot Number: 0.0022 
| | btu ` 
Average Uncertainty in h: IO == 
pr EF 
Time (Secönds) Temperature (?F) h DEU _ 
<a ——— — 9 — — Ee IE 
0.0 ILS) ate 
O NI 38.06 
200 IPPO 3999 
15.0 105.7 40.30 
20.0 JO H 40.32 
2D EE 40.70 
30.0 95.4 40.57 
33.0 91.4 43.12 
575 91.0 2103 
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TABLE VI 


RESULTS OF THERMOCOUPLE CALCULATIONS, HIGH TUNNEL SPEED 


Run: #8 
Tunnel Air Velocity: 306.0 —t 
Sec 
Average Biot Number: 0.0022 
s ; btu 
Average Uncertainty in h: 201296 m — 
Dre E 
Time (Seconds) Temperature (°F) ee 
—— a w EE ss — ii HEER 
0.0 11970 
S 114.4 | 9957 
10.0 LLORO 239.39 
VERD I 105.8 40.81 
20 9 102.2 40.61 
2500 98.6 Zu ds 
3080 95.5 41.23 
35 0 | 32.5 41.70 
38.0 9 1: O 41.53 
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TABLE 12 


RESULTS OF INFRARED SYSTEM CALCULATIONS, LOW TUNNEL SPEED 


Fun: #1 
Tunnel Air Velocity: 201.3 Et 
l sec 
Average Biot Number: 0.00164 
Average Üncertainty in h: = 1.88 — 
hr-ft -°F 
Time (Seconds) Temperature (°F) h (no a 
dee — hr-ft-® 
oe ` la 26.5 
6. 1 14255 3326 
EO 1075 332216 
22.0 102 25 21 29 
yı 5 97. 5 30790 
41.7 ° 2 31 315 
52.5 975 31.99 


45 


TABLES 


RESULTS OF INFRARED SYSTEM CALCULATIONS, LOW TUNNEL SPEED 


Run: #2 
Tunnel Air Velocity: 193 a 
` sec 
Average Biot Number: 0500151 
| ei 
Average Uncertainty in h: t72 EG mo: 7 
I OE 
: weder ' gan — 
Time (Seconds) Temperature (?F) h (Ee) 
Qu LS 
SE 112.5 | 28.40 
bo PEORES 28.46 
os 108,5 27.96 
13% 05 Bs 
1720 102, 5 o 
Di 5 105 272899 
24.4 10025 ZU 
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TABLE 14 


RESULTS OF INFRARED SYSTEM CALCULATIONS, LOW TUNNEL SPEED 


Run: #3 
junnel Air Velocity: More = 
sec 
Average Biot Number: 0 00161 
. : DEU 
Average Uncertainty in h: c 31.2918 5 
ee EF 
Time (Seconds) Temperature (°F) FØR otu ) 
A SEER ca nn a hr-ft-°F 
0.0 E Amo 
5.5 1215230 30.68 
Telks:5 1221250 30852 
19779 10750 29202 
26.5 103.0 28.91 
223 29-0 eS 
40.4 0 31.54 
50.4 91.0 318 32 
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TABLE 15 


RESULTS OF INFRARED SYSTEM CALCULATIONS, LOW TUNNEL SPEED 


Run: #4 
Tunnel Air Velocity: 207.9 Et 
sec 
Average Biot Number: 0.0018 
Average Uncertaint Rina o —2 
emer b 
Time (Seconds) Temperature (°F) h ( btu ) 
= — bii E 
0.0 111 0 
ER 1050 34.61 
6.0 107 0 p 
9.4 10580 32235 
12.4 10580 33.48 
6 103950 32.64 
1903 JON 33.94 
2288 9720 34.44 
290 950 34719 
31.8 93.0 39.65 
56 0 SO 535g 
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TABLE 16 


RESULTS OF INFRARED SYSTEM CALCULATIONS, HIGH TUNNEL SPEED 


Run: #5 
Tunnel Air Velocity: 308.2 —+t 
sec 
Average Biot Number: 099021 
Average Uncertainty in h: nel. E 
r= f reo pP 
Time (Seconds) Temperature (°F) E ( in ) 
+ rn hr-ft-°F 
EO "119.9 
4.7 1150 35.91: 
8.9 11000 39.44 
liði? KOTO 40.08 
193 LOBO 39.70 
29 99.0 42. 7 
296 9500 43.05 
3728 91.0 41.75 
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TABLE l! 7 


RESULTS OF INFRARED SYSTEM CALCULATIONS, HIGH TUNNEL SPEED 


Run: #6 
; ft 
Tunnel Air Velocity: 303.8 — 
sec 
Average Biot Number: 0.0018 
Average Uncertainty in h: BOUES — 
ME — b 
Time (Seconds) Temperature (°F) h T LN 
—  — — AX AE Bt-?F 
0.0 | 119.0 
539 Lise 35.45 
10.8 111 0 3769 
1682 107.0 sc ol 
21005 109590 35.47 
26 oord 37. 
32:19 Eie 38 9 
39.6 Bed 39.86 
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TABLE 18 


EL UCTS OE TNEEARED SYSTEM CALCULATIONS, HIGH TUNNEL SPEED 


Bean: #7 
Rommel Air Velocity: 306.0 ZE 
sec 
Average Biot Number: 0. 0021 
Average Uncertainty in h: mel 06 A 
HE EF 
Time (Seconds) Temperature (°F) h a er 
ar e ee hr-ft-°F 
0.0 11.9559 
5.0 115.0 | 36.14 
9... jp) 38.26 
14.0 INA. > 38.45 
18.9 "93.0 39.70 
24.0 99.0 41.00 
29.7 SL) 41.90 
3658 91.0 41.81 
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TABLE 19 


RESULTS OF INFRARED SYSTEM CALCULAFIONS, HIGH TUNNEL SPEED 


Run: #8 
Tunnel Air Velocity: Fn 
sec 
Average Biot Number: 00022 
Itu 
Average Uncertainty in h: «85 REECH 
hr-ft -?F 
Time (Seconds) Temperature (?F) h as 
— -—— EER nere Et 
049 ee) 
4.9 liao" ` 37.14 
10483 Jelang 37.46 
144% O 5 25 
20.6 103.0 | 37.20 
2579 IO 59735 
Semet» DD 40.45 
30 91.0 41. 32 
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TABLE 20 


RESULTS OF TVC VANE CALCULATIONS 


Run: #1 
Tunnel’ Air Velocity: 2143 E 
Sec 
Average Biot Number: 0.0041 
Average Uncertainty in h: in 3223 a 
hr-ft -°F 
Time (Seconds) ener vet (AF) h ( EL ) 
Pe —— — ——— hr-ft-°F 
0.0 S9 
Deb 95.0 33552 
1.08 91.0 41.41 
1280 S 41.89 
20% 83-0 43.19 
316 79.0 143331 
42.2 75-0 42.50 
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RESULTS OF TVC VANE CALCULATIONS 


Run: 
Tunnel Air Velocity: 
Average Biot Number: 


Average Uncertainty in h: 


Time (Seconds) 


16.4 
2322 
339 


AS 


TABLE 21 
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89 


85 
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#2 


22015 


0590953 


t 4.41 


Temperature. (BN) 


0 


„0 


. 0 


zu 


0 


ji 


0 


54 


Sec 


btu 


hr-ft?-op 
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TABEE 22 


RESULTS OF TVC VANE CALCULATIONS 





Run: #3 
i : eh 
Tunnel ASPEVOTOGCHS 220.5 —— 
sec 
Average Biot Number: 0.0022 
ie 
Average Uncertainty in h: > VU m > 
: I Me Et -°E 
Time (Seconds) Temperature (°F) ner zu 
AAA en ———— iS 
0.0% 97.20 
6.2 50 41.43 
13.0 89 50 42.25 
Zar, Go O 40.95 
Do. 7 81.0 | 45.06 
40.4 IO 44.39 
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APPENDIX D 
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Figure 1. Diagram of Wedge Model 
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Rear View of Wedge Model Assembly 
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Figure 
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Figure 7. Color Thermogram of Wedge 
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Photograph of TVC Vane 


Figure 8. 


Flow Direction 
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Figure 9. Color Thermogram of TVC Vane 
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